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Abstract 

It  haa  recently  been  shown  that  the  physical  properties  of  a  layered 
ocean  floor  may  in  principle  be  determined  from  an  analysis  of  the  reso¬ 
nances  In  frequency,  or  in  the  angle  of  Incidence,  which  appear  in  bottom- 
reflected  acoustic  signals.  In  the  present  paper  we  demonstrate  how, 
using  sinusoidal  signals  of  long  but  finite  duration,  individual  reso¬ 
nances  can  be  selectively  excited  and  analyzed.  Multiple-bounce  layer 
resonances  interfere  destructively  with  the  specular  echo,  and  manifest 
themselves  through  a  characteristic  ringing. 


1  Part  of  this  work  was  carried  out  while  the  authors  were  visitors  to 
the  College  of  Marine  Studies,  University  of  Delaware,  DE  19958. 
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It  has  recently  been  shown  that  the  physical  properties  of  a  layered  ocean  flooi 
nay  in  principle  be  determined  from  an  analysis  of  the  resonances  in  frequency, 
or  in  the  angle  of  incidence,  which  appear  in  bottom-reflected  acoustic  signals. 
In  the  present  paper  we  demonstrate  how,  using  sinuosoidal  signals  of  long  but 
finite  duration,  individual  resonances  can  be  selectively  excited  and  analyzed. 
Kultiple-bounce  layer  resonances  interfere  destructive  with  the  specular  echo, 
and  manifest  themselves  through  a  characteristic  ringing. ^ 
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1.  Introduce Ion 

Remote  acoustic  sensing  of  the  ocean  floor  conventionally  employs 
pulse  arrival  time  methods ^see,  e.g.,  Houtz,  1974;  Bell  and  Porter, 

1974;  Maynard  et  al.,  1974^,  which  permit  a  determination  of  the  ocean 
floor  depth  profile  and,  to  some  extent,  also  of  the  thicknesses  and 
sound  velocities  of  sediment  layers  but  not  of  the  layer  densities.  We 
recently  proposed  a  new  method  for  the  remote  sensing  of  ocean  floor 
properties  ^Nagl,  Uberall  and  Hoover,  1982  j which,  at  least  in  principle, 
may  lead  to  a  determination  of  all  layer  parameters  by  remote  acoustic 
sensing,  provided  the  layer  structure  is  well  defined.  This  method  is 
based  on  analyzing  the  effects  of  layer  resonances  which  are  carried  by 
acoustic  signals  reflected  from  the  ocean  floor;  and  as  shown  by  us.  It 
relates  measured  resonances  parameters  (locations,  heights,  widths) 
explicitly  to  the  medium  parameters  (densities,  thicknesses,  sound  veloci¬ 
ties  etc.  of  the  layers),  hence  providing  a  solution  to  the  inverse 


scattering  problem. 

In  the  mentioned  reference  ^Nagl,  Uberall,  and  Hoover,  1982^,  the 
steady-state  reflections  of  purely  sinuosoldal  signals  were  considered,  as 
applied  to  the  case  of  a  fluid  layer  sandwiched  between  two  different  fluid 
half-spaces  with  the  impedances  increasing  with  depth,  this  case  being 
the  simplest  model  relevant  to  a  sediment-covered  ocean  floor.  Using  the 
same  model,  we  here  study  the  implications  of  employing  finite-length 
sinusoidal  wave  trains  as  acoustic  probes.  These  enable  us  to  directly 
study  individual  layer  resonances  and  to  observe  their  ringing,  thus 
directly  measuring  the  resonances  width.  Furthermore,  they  permit  a  separa¬ 
tion  of  the  reflected  signal  into  a  specularly  reflected  component,  and  one 


\ 
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that  has  undergone  multiple  reflections  within  the  layer,  as  well  as  a 
determination  of  the  interference  between  the  two  components.  Consequently, 
more  information  is  obtained  from  the  use  of  finite-length  wave  trains 
rather  than  of  steady-state  sinuosoids.  In  fact,  it  will  be  shown  that 
in  principle  two  long  wave-trains  reflected  at  different  angles  contain 
enough  information  to  permit  a  determination  of  all  five  bottom  parameters 
in  the  three-fluid  model  considered  here. 


2.  The  Reflection  Coefficient  for  the  Three-Fluid  Model 


The  geometry  of  the  three-fluid  model  is  shown  in  Fig.  1.  Sound  is 
incident  through  the  ocean  water  (mediun  3),  reflects  from  or  penetrates 
the  sediment  (medium  2),  then  reflects  from  or  penetrates  the  substratum 
(medium  1).  The  sound  speed,  density,  and  propagation  direction  for  the 
ith  medium  are  designated  c^,  p^,  and  0i,  respectively,  and  d  is  the 
thickness  of  the  sediment  layer. 

The  total  reflection  coefficient,  as  a  function  of  wave  number  k”w/c^, 

Nagl,  Uberall  and 

Hoover,  1982 J 


*oo 


a  cos  r  -  4  (i  -  £ 

Sf  cos  -  <  ( 1  *  r1') 


(la) 


where 


*cLto 


cos  9X 


X  *■- 


(lb) 


containing  the  impedance 


to /co*0*  <ld 

of  the  ith  medium.  This  reflection  coefficient  may  be  decomposed  into 
the  form 


where 


and  where 


^  cot  e,  , 


2-?, 
Z,*2  ’ 


£  » 1,Z 


are  the  reflection  coefficients  of  waves  propagating  in  layer  2  and  being 
reflected  from  the  boundaries  with  the  media  1  and  3,  respectively.  Geo¬ 
metrically,  Eq.  (2a)  corresponds  to  the  multiple- ref lection  picture  shown 
in  Fig.  2.  The  phase  factors  exp  ikOin  Eq.  (2a)  account  for  the  phase 
shift  Introduced  by  one  round-trip  bounce  in  layer  2.  Sunning  Eq.  (2a) 


leads  to 


m  = 


>  kt> 


«  t  -  A( s>V|,v nr  (2d) 

fa+t)wi  (l  +  *i>)  SiV 

with  a  ■  Z-JZjt  b-zyZj,  equivalent  to  Eq.  (la). 

In  Fig.  3,  we  plot  J  R(lt) | ^  vs.  the  dimensionless  frequency-thickness 
product  fd  (where  f»U>/2‘7r )  at  a  fixed  angle  of  incidence  0^  (here  taken 
as  normal  Incidence,  9^  “  0°),  for  the  model  parameters 

P-  ■  6.5,  P-  •  3.3,  P  -  ■  1.0  g/cm'* 

^1  2  3  (3) 

cx  -  5495,  c2  -  2544,  c3  -  1500  m/s. 


This  graph  shows  "antiresonances"  spaced  at  regular  frequency  Intervals. 

A  similar  picture  is  obtained  if  |r|*  is  plotted  vs.  at  fixed  k, 
except  that  here  the  (anti) resonances  have  different  width  and  are  non- 
uniformly  spaced.  The  layer  parameters  were  chosen  such  as  to  achieve 
high  impedance  contrasts  at  both  interfaces  in  order  to  emphasize  the 
resonance  characteristics  of  the  reflected  signal.  If  smaller  (more 
commonly  encountered)  bottom  parameters  were  used,  the  resonances  in 
Fig.  3  could  appear  less  pronounced,  but  all  results  obtained  below  still 
apply. 

The  resonance  condition  corresponds  to  the  vanishing  of  the  denominator 
of  Eq.  (2d) 

satisfied  by  the  complex  resonance  wave  numbers 
*  (2m*  1 )ie  -  < 

CD  *  ^  (l/PiR,), 

whose  real  parts  determine  the  resonance  frequencies 

«  “Cj  Re  kw  *  'Cj  • 

At  the  resonances  one  thus  finds 

4+f  4  C  ^  ^  * 

the  sign  being  determined  by  the  requirement 

«*P 

Since  in  the  ocean  floor  the  impedances  generally  increase  with  depth 
(Z^>  Z2^Z^),  one  has  Rj>  0,  R^<  0  and  the  resonances  are  found  at 

le^T>  »  *  (2*i+l)  ir,  , 


(4) 

(5a) 

7  0, 

(5b) 

• 

• 

n 

(5c) 

(6a) 

(6b) 

l.e.  for 


2-w*  -d 

2d.  0^ 


(7a) 


(7b) 
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or 


4^  - 


(7c) 


Applying  Snell's  law  puce  Che  resonances  ac  Che  values 
D  j  _  (2-w-rl)^ 

of  Che  frequency-thickness  produce. 


(8a) 


From  Eq.  (6a)  lc  Is  now  apparenc  chac  ac  resonance  all  multiple- 
bounce  eerms  in  R(k),  Eq.  (2a),  add  coherenely.  Remembering  chat 


(8b) 


R^<0,  we  obtain  che  reflection  coefficient  at  resonance: 

fcflo-  iw-  **  nw  • 

Here,  a  large  term  is  subtracted  from  che  specular-reflection  component 
|  R^j ,  leading  to  the  minima  ("anti-resonances")  J R^Ck^) J  2  in  Fig.  3. 

In  contrast,  the  denominator  of  Eq.  (2d)  reaches  a  maximum  for 
exp  ikD  »1.  The  associated  frequencies 

&xL  “  <>c> 
lie  halfway  in  between  the  resonance  frequencies.  Successive  multiple 

reflections  add  destructively  (with  alternating  signs)  in  Eq.  (2a),  and 
in  che  corresponding  reflection  coefficient 

i t4S:\  <M) 

a  small  term  is  subtracted  from  |  | ,  leading  to  the  broad  maxima  in  Fig.  3. 

The  resonances  at  f  are  thus  transmission  resonances,  or  anti-resonances  in 
n 

the  reflection  coefficient,  physically  characterized  by  che  multiple  reflec¬ 
tions  adding  in  phase,  but  their  sum  interfering  destructively  with  the 
specular  reflection  amplitude. 
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3.  The  Reflected  Signal  for  Incident  Sinusoidal  Have  Trains 


Finite-length  wave  trains  nay  be  used  to  separate  specularly  and 
aultiply-reflected  signals,  due  to  the  differences  in  arrival  tines 
Involved.  A  sinusoidal  wave  train  of  duration  tp  -  2irnQ/u)0  (nQ  being 
the  number  of  cyles  in  the  pulse  which  we  shall  assuae  integer)  and  carrier 
frequency  U).  *  kQ  has  the  spectrum 

e  *  *  -  i 


G(k)=  k. 


(9) 


k%-  ko*- 

peaked  around  kg  or  fg  within  the  intervals  of  wave  number  k  or  frequency  f: 

A K-  2k./-n.  ,  Zf./-n.  »  (10a) 

respectively.  If  now  this  interval  is  small  compared  to  the  resonance 
spacing 


i .  e. , 


^  ^  2a.[('C,/c4)t-si'vl,0>T/*  ’  (1°b) 


z£»/'w.#  46-  /C z/2  d  cos9x  , 


(lla) 


and  if  one  lets  the  carrier  frequency  f^  coincide  with  a  resonance  frequency. 


i.  -  u  - 


-  2.»\+i 


4*®,©^  S  any 

one  will  selectively  excite  this  Individual  resonance.  Equations  (lla) 


and  (lib)  give 

nQ^2n  +  l,  (He) 

determining  the  required  minimal  length  of  the  incident  pulse 

<»,*,* )*‘fC(k)eik("h‘*>  W*c.  (12.> 

—  »o 

Using  for  the  reflection  coefficient  in  the  reflected  pulse  , 

p„f  (*,  »,t)  -  (ecu)  j  y/*r,  (12b) 


.  1 
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(13a) 


(13b) 


the  expression  of  Gq.  (2a)  rewritten  as 

1200  =  £  E* 

4MSO 

Xo  *  "  » 

J>J.* 

XL*  ** 

Cell lx^ the  observer's  position  (Xg,  zQ)  and  using  a  time  variable  t^  which 
vanishes  when  the  leading  edge  of  the  reflected  pulse  reaches  the  observer, 
one  finds 

iw  *. ,-o  -  E  a,  \  ©  (*,  -  if ) «« -  <**  -  S* ) 

'**  (14a) 

-  ©  fa  ~  ^  5iV  u>0  ( t  j  ^  j  . 


where 


@(''r) 


ro 

*(1  r> 


TT<-  O 
O 


(14b) 


is  the  Heaviside  function.  This  corresponds  to  a  series  of  overlapping 
wove  trains  whose  outline  is  depicted  in  Fig.  4.  The  first  rectangle  outlines 
the  specularly  reflected  pulse,  the  following  ones  the  pulses  that  underwent 
(multiple)  internal  reflections.  Since  their  coefficients  and  have 

equal  signs  for  m**0  and  opposite  signs  for  ra>  0,  and  since  at  every  resonance 
the  relative  phase  difference  between  any  two  successive  sine  functions  in 
Eq.  (14a)  can  be  shown  to  be  1C,  one  sees  that  the  individual  wave  trains 
are  in  phase  with  each  other  for  m  ? 0,  but  are  out  of  phase  with  the 
specular  pulse  (art)). 

For  the  example  with  layer  parameters  given  by  Gq.  (3),  and  assuming 

normal  incidence  throughout,  we  evaluate  numerically  the  reflected  signal 

of  Iq.  (14a).  The  resonances  spacing  is  Af  ■  (1272/d)He  ,  with  d  measured 

in  meters.  We  choose  a  spectral  width  AFa0.1Af  so  that  Hq  ■  20fg/A  f 

and  t  •  20/ A f,  and  let  the  carrier  frequency  coincide  with  the  third 
P 


Vs* 


J 


resonance  (n*2)  so  that  fo  *  2.5Af.  Figure  5  shows  the  reflected  signal, 
starting  at  the  receiver's  time  With  the  Internally  reflected 

components  coming  in  at  intervals  of  D/c^  (here  equal  to  2.5  cycles  of  the 
carrier),  the  signal  is  steprwise  reduced  to  its  steady-state  amplitude. 

After  the  passing  of  the  trailing  edge  of  the  specular  pulse  at  t^t^,  only 
the  internally  reflected  pulses  remain,  manifesting  themselves  as  a  ringing 
of  the  selectively  excited  resonance. 

If  the  carrier  frequency  is  chosen  halfway  between  two  adjacent  reso¬ 
nances,  the  results  are  as  shown  in  Fig.  6  for  a  choice  of  fQ*3.0Af. 

Successive  internal  reflections  interfere  destructively,  causing  the 
alternating  steps  in  the  initial  transients  leading  to  the  steady-state 
amplitude.  Some  ringing  is  seen  to  persist  even  here  (at  a  small  amplitude 
level),  since  finite  resonance  widths  cause  some  overlap  in  the  resonance 
peaks. 

Figure  7  shows  the  case  of  fQ  »  2.6£f,  i.e.  the  carrier  frequency  \ 

being  slightly  off-resonance  (by  1/10  of  the  resonance  spacing).  The  steady- 
state  amplitude  has  already  changed  significantly  compared  to  the  on-resonance 
case  of  Fig.  5,  indicating  a  narrow  resonance  width. 

The  half-width  of  the  resonance  peaks  (which  is  here  the  same  for  all 
resonances,  cf.  Fig.  3)  is  related  to  the  imaginary  part  of  k  ,  Eqs.  (5),  by 


(case  RjRj<0).  Since  fQ  "  (n+*s)Af,  we  have 

Y*  _  Af  D  (  d  ^ 

T  -  >  ■  '15k> 

The  ringing,  seen  from  Eq.  (14a)  and  Fig.  5  to  take  place  in  a  step-wise 


decaying  fashion  at  the  time  Intervals  A  t«l/A  f  -  D/c^: 


,  ^ .  j-—  ’  4>V 


(16a) 


J 
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may  be  approximated  by  a  decaying  exponential 

P  ref(t-tp4nAt)  »  t>tp*  (16b) 

Comparison  gives 

^  -Af A  (-1/R1R3) ,  (16c) 

or 

Jr  r  '  (16d) 

which  shows  that  the  resonance  half-width  may  be  directly  determined 

by  measuring  the  decay  content  of  the  ringing.  From  Eq.  (8a)  it  is 
apparent  that  the  resonance  spacing  is  larger  for  oblique  incidence  than 
from  normal  incidence  (provided  that  ,  which  is  assumed) .  The 

resonance  width  as  given  by  Eq.  (15b)  may,  on  the  other  hand,  increase  or 
decrease  with  the  angle  of  incidence.  However,  the  relative  halfwidth 
(the  halfwldth  divided  by  the  resonance  spacing)  always  decreases,  i.e. 
the  resonances  appear  more  pronounced  for  oblique  incidence  than  for 
normal  incidence. 


4.  Information  Content  of  the  Reflected  Signal 

For  a  steady-state  incident  signal,  the  only  information  obtained  is 
that  contained  in  the  total  reflection  coefficient. 


(17) 


With  a  finite-length  sinusoidal  pulse  train,  however,  one  may  measure 
separately  the  three  quantities: 

(the  amplitude  of  the  specular  reflection). 


! 
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1  (i  +  R^Rj1  -  2^^  C*. 


(18b) 


or,  more  conveniently. 


a1a>  -  Bj(i  -  Rj2, ) 


(18b’) 


(the  amplitude  of  all  internal  reflections,  or  only  of  the  first  one, 
respectively),  ajJ 

At  -  1/Af  -  Cj/D  (18c) 

(the  delay  time  between  successive  internal  reflection).  The  measurement 
of  Rj  or  Rj^  can  be  seen  to  be  equivalent  to  that  of  ^  ,  cf.  Eqs.  (16). 
From  our  previous  analytic  expressions  of  these  quantities,  one  may  now 
derive  three  expressions  that  relate  Rg,  and  to  the  layer  parameters: 


and 


jL^  _  i~%S  3fj  A  , 
i+Ki  a  M> 

i  _  (  A>t  \x  (  riV0>  Y" 

'Cz  \  "l  ^  1  > 

-  ft  -  (*»*?  . 


(19a) 


(19b) 


(19c) 


where  we  defined 
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i-x:  -Rj  i-Kj  • 


(19d) 
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This  constitutes  three  equations,  Eqs.  (19a)-(19c),  for  the  five 
unknowns  cj_»  c2»  fi  $2  and  d  °f  our  ocean  floor  model.  Information 
on  the  resonance  order  n.  Incidentally,  is  obtained  from  Eq.  (lib), 
since  the  resonance  spacing  A  f  is  known  from  the  measurement  of  At  , 
Eq.  (18c). 

Two  Independent  measurements  of  Rs,  and  At  at  the  same 
frequency,  but  at  two  different  angles  of  incidence  and  63^ »  can 

now  be  shown  to  lead  to  sufficient  additional  information  so  that  one  may- 
solve  the  resulting  equations  for  all  bottom  parameters  involved.  Indeed, 
Eq.  (19b)  then  yields  directly 


*>[  aC-  A-t,1  y 

■2  U-* 

f  h-ti-M? 

‘  °3  (' 


r. 


(20a) 


(20b) 


furnishing  the  unknowns  d  and  Cg.  With  d  known,  Eq.  (19a)  can  be  used  to 
find  ^>2’ 

Similarly,  Eq.  (19c)  yields  the  remaining  two  unknowns, 

*  _  )  > 


( 


V 


<»> 


(20c) 


fi=  [  (2T/01  -  &,*  Y\ 


t  (20d) 


In  the  above  equations  ,  A t^,  8^3^^  and  are  the  values  of  At, 

and  Zj  pertaining  to  the  angles  9^,  respectively.  Thus,  when  long  but 
finite  sinusoidal  sound  pulses  are  used  as  probes,  in  principle  only  two 


J 


-In¬ 


different  measurements  of  Rg,  (or^T  )  and  At  have  to  be  made  at  the 
same  frequency  and  at  tvo  different  angles,  in  order  to  deduce  all 
unknown  layer  parameters. 
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The  acoustic  resonance  method  for  the  remote  sensing  of  ocean  floor 
properties  was  introduced  in  our  previous  work  ^Nagl,  Uberall  and  Hoover, 
1982 J  for  the  case  of  monochromatic,  steady-state  sinusoidal  signals  used 
as  probes.  In  that  case,  resonances  of  ocean  floor  sediment  layers,  e.g., 
in  frequency  (or,  alternately,  in  the  angle  of  Incidence)  must  be  found  by 
repeating  the  experiment  at  many  different  values  of  the  frequency  (or  of 
the  angle).  With  the  present  suggestion  of  using  finite-length  sinusoidal 
wave  trains,  the  number  of  measurements  can  be  dramatically  reduced.  The 
pulse  length  must  be  sufficiently  great  so  that  the  specularly  reflected 
amplitude  and  several  of  the  signals  multiply  reflected  from  within  a 
bottom  layer  overlap;  it  is  then  possible  to  selectively  excite  and  study 
one  individual  resonance.  The  time  delay  between  internal  reflections 
determines  the  order  of  the  resonance  and  the  resonance  spacing.  Measure¬ 
ment  of  specularly  and  internally  reflected  amplitudes,  or  of  the  ringing 
of  the  resonance,  carried  out  at  tvo  different  angles  then  suffices  in 
principle  to  furnish  us  with  the  values  of  the  layer  and  substratum 
parameters,  hence  solving  the  Inverse  scattering  problem.  While  all  this 
was  demonstrated  here  on  the  simple  example  of  a  fluid  layer  on  top  of  a 
fluid  substratum,  it  is  to  be  expected  that  a  multilayered  ocean  floor  will 
in  principle  yield  similar  information  due  to  layer  resonances  beitg  set  up 


in  individual  layers  by  an  incident  finite-length  wave  train.  Work  on 
this  problem  is  now  in  progress. 

It  should  be  noted  that  observations  of  the  ringing  of  resonances, 
acoustically  excited  by  finite-length  wave  trains,  have  been  used  [Maze, 
Taconet  and  Ripoche,  198  l^J  in  order  to  determine  the  resonance  spectrum 
and  the  resonance  widths  of  metal  cylinders  iamersed  in  a  fluid.  These 
observations  have  been  interpreted  [Mumrich  et  al.,  1984 3  in  terms  of 
overlapping  wave  trains  circumferentially  propagating  around  the  target, 
in  a  fashion  similar  to  the  multiply-reflected  wave  trains  of  the  present 
paper.  Likewise,  overlapping  wave  trains  in  radar  scattering  have  been 
analyzed  analogously  [towell  and  Uberall,  1984 
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Fig.  2 

Fig.  3 

Fig.  4 
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.  Geometry  of  tbs  three-fluid  model  for  an  ocean  floor  with  a 
sediment  layer. 

.  Geometrical  view  of  multiple-reflection  decomposition  for  the 
reflection  coefficient  of  a  sediment  layer. 

.  Plot  of  | R(k)|  ^  vs.  frequency-thickness  fd,  at  normal  incidence 
(©3  ■  0°),  showing  (anti) resonances. 

.  Amplitudes  of  reflected  overlapping  wave  trains. 

.  Reflected  signal  vs.  time,  for  the  carrier  frequency  coinciding 
with  third  resonance  frequency.  Layer  parameters  of  Eq.  (3). 

.  Reflected  signal,  for  the  carrier  frequency  halfway  between  the 
third  and  the  fourth  resonance  frequency. 

Reflected  signal,  for  the  carrier  frequency  lying  above  the 
third  resonance  frequency  by  1/10  of  the  resonance  spacing. 
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